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Abstract. - The phenomenon of intrinsic localization in discrete nonlinear extended systems, 
i.e. the (generic) existence of discrete breathers, is shown to be not restricted to periodic 
solutions but it also extends to more complex (chaotic) dynamical behaviour. We illustrate 
this with two different forced and damped systems exhibiting this type of solutions: In an 
anisotropic Josephson junction ladder, we obtain intrinsically localized chaotic solutions by 
following periodic rotobreather solutions through a cascade of period-doubling bifurcations. 
In an array of forced and damped van der Pol oscillators, they are obtained by numerical 
continuation (path-following) methods from the uncoupled limit, where its existence is trivially 
ascertained, following the ideas of the anticontinuum limit. 



Discrete homogeneous arrays of (hamiltonian and non-hamiltonian) nonlinear oscillators (or 
rotors) exhibit generic solutions which are time-periodic and (typically exponentially) localized 
in space. These solutions are called discrete breathers by analogy with non-topological local- 
ized solutions of certain PDE's. In contrast with continuous "bona fide" breathers, discrete 
breathers posses a remarkable structural stability, and thus genericity. This localization is 
often referred to as intrinsic to stress the fact that the system is homogeneous (no impurities 
or disorder are present). For an updated and comprehensive review on discrete breathers, 
see|. 

A general schematic way to describe a discrete breather in a one-dimensional lattice is the 
following: Let us consider the phase space T s of a single oscillator (or rotor), so that the phase 
space T of the network is the cartesian product of the single site phase spaces. Let denote 
by A, B periodic orbits in T s , eventually projections of trajectories of Y onto T s . A discrete 
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breather is a solution 

{4> i {t)} = {...,B_ 2 ,B_ 1 ,A,B 1 ,B 2 ,...} (1) 

with limiji^oo Bi = B^, and B^ ^ A. Archetypical examples are Klein-Gordon hamiltonian 
breathers, where A is a periodic cycle of frequency Wf, in the (4>, <p) phase space, B x is the rest 
solution (0,0), and Bi are Wb-cycles with exponentially decreasing amplitude. In the case of 
forced and damped arrays, A and B^ are usually u^-cycles of different amplitude. If A is an 
Wb-cycle non-homotopic to zero (i. e. the central oscillator rotates), the term rotobreather is 
used§,|. 

In this letter we present numerical evidences as well as plausibility arguments strongly 
supporting the conclusion that the phenomenon of intrinsic localization in discrete nonlinear 
extended systems is not restricted to time-periodic solutions, but it extends to more complex 
(chaotic) behaviour in a generic way for damped and forced systems. More specifically, we show 
below examples of solutions of the type schematized in (Q), where A is a chaotic trajectory, 
Boo is a "regular" a^-cycle, and Bi are "noisy" cycles, with "noise intensity" exponentially 
decreasing to zero as \i\ grows. The first example concerns the operation of a Josephson 
junction device, the Josephson junction ladder, which has recently received some attention 
from both theoretical j|, [j| ^, 0| and experimental || sides, in connection to the relevance of 
nonlinear dynamics of discrete systems in Condensed Matter Physics. The second example, 
though also experimentally realizable, serves us to illuminate possible pathways towards a 
rigorous characterization of the gencricity of intrinsically localized chaos in discrete nonlinear 
extended systems, in the spirit of the ideas of the, so called, anticontinuum limim\ IToL 11 



approach to intrinsic localization. We end with a short discussion on the implausibility of 
existence of this type of solutions as exact ones in hamiltonian arrays. Earlier numerical 
observations of localized chaotic solutions seems to have been reported in|12] (coupled map 



lattices) and|13 (domain walls in a parametrically excited lattice of oscillators). Our results 
establish a precise (and very general) link between situations of spatio-temporal complex 
behaviour in spatially extended discrete systems and the emergent new results and powerful 
methods of intrinsic localization. 

Recent theoretical analyses of the dynamics of an anisotropic Josephson junction ladder 
(see figure [l]) with injected ac currents (5) have shown the existence of discrete breathers as 
attracting solutions of the equations of motion describing the dynamics of the system in the 
framework of the resistively and capacitively shunted junction (RCSJ) approach |l4|]. The 
existence of discrete breathers in Josephson junction arrays should indeed be regarded as 
generic, given the connection between the general description of these systems in terms of 
the superconducting Ginzburg-Landau order parameter ty(x) = \ty(x)\exp(iO(x)), where x 
denotes the superconducting island position, and the discrete nonlinear Schrodinger equation, 
for the case of ideal (perfect insulating) junctions p]. In fact, the quantum Hamiltonian 
of a single ideal Josephson junction corresponds to the problem of two coupled anharmonic 
quantum oscillators, for which the asymmetric classical breather solutions have been shown to 
persist in the quantum regime as very long lifetime states p5[ (see also Jl6{). When the energy 
cost to add an extra Cooper pair on a neutral superconducting island {charging energy E c ) 
is much lower than the tunneling energy (Josephson energy Ej) the superconducting phase 
9(x) becomes a good (very weakly fluctuating) variable to describing the island state, thus 
validating the RCSJ approach JlJ] . This is the situation when the superconducting islands are 
of macroscopic size. The validity of the RCSJ approach in the regime E c /Ej C 1 is a well 
established issue and its predictions fit very well with experiments |18[ . 

0i and 0'i will denote, respectively, the phases of upper and lower islands at site i in the 
ladder; the currents I(t) — I ac cos(wi) are injected into the islands in the upper row and 
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Here is the figure 



Fig. 1. - Schematic picture of the J J ladder showing the injection of the currents in the array. 

extracted from those in the lower row; (J x , e x ) are the junction characteristics for junctions 
in horizontal links and (J y , e y ) for junctions in vertical links. With the change of variables 
Xl = I(6». ( + O'i), fa = ±(0* - 6K), the RCSJ equations [§ are 

Xi = J x [sin(xi+i - Xi) cos(fa +1 - fa)+ sin(xi-i - Xi) cos(fa-! - fa)] 

+e x (Xi+i + Xi-i - 2%i) (2) 

fa = Jx [cos(xi+i - Xi) sm(fa+i - fa) + cos(x 4 -i - Xi) sin(fa-i - fa)] 

+e x (j>i+i + fa-i - Ifa^j - J y sin(2^) - 2e y fa - I(t) (3) 

With uniform initial conditions in the "center of mass" coordinates and momenta: Xi an d 
Xi independent of i, equations (^) have the solution Xi(i) = fit + a for all i; this effectively 
decouples equations (|^) for the fa variables from equations (|^) for the Xi variables. Then, using 
efficient continuation methods [ [To| , pj| from the uncoupled (anticontinuum) limit ( J x — e x =0), 
one easily computes discrete breather solutions; these turn out to be attractors of the dynamics 
of the ladder in a wide range of parameter values. 

We will concentrate on the rotobreather type of solutions, in which the phase half-difference 
<pj* through a vertical junction rotates, while the rest fa (i ^ j*) oscillate, and the "center 
of mass" variables Xi remain uniformly at rest (f2 = a = 0; note that any other values for 
these parameters, fixed by the uniform initial conditions, would show the same behavior). The 
period of the rotobreather solution is = 2n/ujb — 47r/w, where uj is the frequency of the 
external currents. 

By performing the Floquet analysis of rotobreather solutions, one can determine the re- 
gions of linear stability in parameter space, whose borders correspond to different types of 
bifurcations 0. One of them (which occurs typically when varying the external frequency 
oj) is a period-doubling bifurcation: The (destabilizing) eigenvector of the Floquet matrix, 
which is associated to the eigenvalue exiting the unit circle (in complex plane) at —1, is 
(exponentially) localized at the center of the rotobreather and then, a new (linearly stable) 
rotobreather with frequency ujb/2 exists past the bifurcation. This new rotobreather can be 
easy and safely obtained by slightly perturbing the unstable rotobreather along the direction 
of the destabilizing eigenvector. In other words, although one cannot continue the localized 
solution in a bifurcation, local bifurcation analysis helps to throw a bridge over the bifurcation, 
so arriving safely to the new localized solution at the other side. 

Continuously varying the external frequency lu, further period doubling bifurcations are 
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Here is the figure 



Fig. 2. - Lyapunov spectrum of a chaotic rotobreather in the JJ ladder. 



often found leading to a chaotic solution. In order to characterize unambiguosly this solution 
as chaotic, we have computed its Lyapunov spectrum {A^}, which is shown in figure |^. There 
is only one positive Lyapunov exponent, Ai = 0.049bits/s. As we are dealing here with a 
continuous time dynamical system, a null exponent is also present. The rest of the spectrum 
is negative. Thus, there is only one expanding direction (degree of freedom) in phase space. 
The estimated Lyapunov dimension, Dl, defined [{17J as 



i=l 



with j such that J2i=i Ai > and J2l=i \ < (exponents are ordered in decreasing order), 
is£> L =4.7. 

A look at the profile (at different times) of the Lyapunov vector associated with the positive 
Lyapunov exponent reveals that it is strongly localized in space. As the period doubling 
bifurcations leading to the chaotic solution are driven by exponentially localized eigenvectors 
of the Floquet matrix, it is not surprising that this chaotic solution is exponentially localized. 
In figure [| we show the Poincare (stroboscopic, with period 2T&) section of the central rotor 
trajectory <po(t) (mod 2tt) of the intrinsically localized chaotic solution for parameter values 
J x = 0.05, J y = 0.5, e x — 0.03, e y = 0.01, w = 1.623 and I ac = 0.72. As shown also in 
figure ||, the trajectories 4>i(t) for \i\ > are noisy (or chaotically perturbed) oscillations. As a 
rough measure of " noise intensity" , we adopt the radius of the smallest circle containing the 
Poincare section of the ith oscillator. This quantity decreases exponentially ~ Cexp(— \i\/£) 
(£ ~ 1.13), as evidenced in figure |[ 

Vaguely speaking, one could say that the uniformly oscillating solution is robust enough to 
exponentially damp out the penetration of the chaotic perturbation produced by the central 
rotor; equivalently, one could say that the uniformly oscillating solution posses finite coherence 
length £, so that an oscillator does not feel the effect of any sustained local perturbation located 
at distances much greater than £ (lattice units) from it. On intuitive basis, it is clear that 



finite coherence length is required for intrinsic localization to occur 1 19 . 

Now we turn to the question on genericity, i. e. should one expect that these intrinsically 
localized chaotic solutions exist generically in discrete arrays of coupled nonlinear oscillators? 
Though arguably there is little doubt that finite coherence length is ubiquitous in discrete 
nonlinear extended systems, at least some degree of robustness of the chaotic trajectory in the 
central oscillator (not to speak of the mere possibility of a chaotic behaviour) is also needed. In 
an attempt to pave the way towards rigorous answers to the question, we have considered the 
perspective on intrinsic localization opened by the " anticontinuum limit" |)], [HJ approach, 



Short Title 



5 



Here is the figure 



Fig. 3. - Poincare (stroboscopic) sections of the trajectories (j>i(t) at times t + n2Tt- The pictures show 
in quite different scales the planes 4> (in radians per unit time) vs. <j> (in radians) for the rotor (a), its 
fifth neighbour (b) and the ninth one (c). In (d), the logarithm of the "noise amplitude" plotted vs. 
the neighboring index shows the exponentially localized character of chaos. 



as explained below. 

Let us consider a chain of forced and damped identical uncoupled oscillators, and assume 
that there is coexistence of a chaotic attractor and an attracting cycle in the single oscillator 
phase space. Now, consider the cartesian product of a (central site) chaotic attractor and 
attracting regular cycles in the rest of lattices sites. This set is an attractor in the phase space 
of the uncoupled chain, which could plausibly be continued when coupling is turned on. 

In order to check this idea, we have chosen a chain of harmonically coupled, forced van der 
Pol oscillators: 

4>i = -fi((t>i - l)4>i - <t>i +bcos(cut) + C{4> i+ x - 2<t>i + <j>i-i) (5) 

For /i = 4.033, b = 9.0 and u> = tt, the single forced van der Pol oscillator phase space shows 
coexistence of two strange attractors and a periodic cycle of frequency ui/3 (see ^(|). We 
have numerically continued the solution of the uncoupled chain in which the central oscillator 
follows a chaotic trajectory in one of the strange attractors, while the rest of the oscillators 
follow the periodic cycle, for non-zero values of the coupling constant C, up to values of the 
order of 0.5 x 10 -3 , which are small but significantly different from zero. 

The continuation from the uncoupled limit provides a systematic way of obtaining intrin- 
sically localized chaotic solutions, provided the coexistence of strange and periodic attractors 
for a single oscillator. It may also serve, like in the simpler case of periodic discrete breathers, 
as a basis for the construction of a proof of existence which we see as a difficult problem. 



Indeed, Mackay |21 already mentioned this approach for the case of the Plykin attractor, 
where continuation is ensured due to uniform hyperbolicity; unfortunately, as usual in chaos 
theory, strong conditions which simplify mathematical proofs do not seem to fit easily into 
realistic physical models. 

The examples we have shown here concern systems of forced and damped oscillators, and 
one may wonder about hamiltonian arrays of oscillators. Though we do not have a definite 
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answer on the existence of intrinsic localized chaos in discrete nonlinear Hamiltonian extended 
systems, it seems plausible that the typical "broad band" structure of the power spectrum 
of chaotic trajectories would imply a violation of the condition of non-resonance with the 
phonons |[o| . In the extent that this condition plays an essential role in the proof of existence 
of hamiltonian discrete breathers, we think that the answer is negative. However, chaotic 
breathers in discrete Hamiltonian arrays easily appear as long-lived transient solutions. An 
observation of erratically moving transient chaotic breathers in hamiltonian Fermi-Pasta-Ulam 
chains has been recently reported (2^]. After completion of this work, we became aware of the 
numerical observation of chaotic rotobreathers by Bonart and Page(2^] in a Id driven damped 
lattice of dipoles. 

We acknowledge to S. Aubry, C. Baesens, R.S. Mackay and J.L. Marin for many useful 
discussions, P. Grassberger for his illuminating criticisms and J. Page for sending us a draft 
of his work, prior to publication. This work has been financially supported by DGES through 
project PB95-0797. One of us (JJM) acknowledges a Fulbright-MEC fellowship. 
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